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We derive an effective Hamiltonian for Gai_iMna,As in the dilute limit, where Gai-ajMn^As can 
be described in terms of spin F = 3/2 polarons hopping between the Mn sites and coupled to the 
local Mn spins. We determine the parameters of our model from microscopic calculations using both 
a variational method and an exact diagonalization within the so-called spherical approximation. Our 
approach treats the extremely large Coulomb interaction in a non-perturbative way, and captures the 
effects of strong spin-orbit coupling and Mn positional disorder. We study the effective Hamiltonian 
in a mean field and variational calculation, including the effects of interactions between the holes 
at both zero and finite temperature. We study the resulting magnetic properties, such as the 
magnetization and spin disorder manifest in the generically non-coUinear magnetic state. We find 
a well formed impurity band fairly well separated from the valence band up to lactive ^ 0.015 for 
which finite size scaling studies of the participation ratios indicate a localization transition, even 
in the presence of strong on-site interactions, where a:activc < a;nom is the fraction of magnetically 
active Mn. We study the localization transition as a function of hole concentration, Mn positional 
disorder, and interaction strength between the holes. 

PACS numbers: 75.30.-m,75.47.-m,75.50.Pp 



I. INTRODUCTION 

Recently there has been a surge of interest in 
the more than 30 year old field of diluted magnetic 
semiconductors^ that has been largely motivated by the 
potential application of these materials in spin-based 
computationSiMik devices. In particular, the discovery of 
ferromagnetism in low-temperature molecular beam epi- 
taxy (MBE) grown Gai_a;Mna;As has generated renewed 
interest.^ In this material Curie temperatures as high as 
Tc w 160K have been observedii 

In this paper we focus on one of the most studied mag- 
netic semiconductors, Gai-jjMn^jAs, though most of our 
calculations carry over to other p-doped III-V magnetic 
semiconductors. In Gai-^jMn^jAs substitutional Mn^"*" 
play a fundamental role: They provide local spin 5 = 5/2 
moments and they dope holes into the latticctS Since the 
Mn^"*" ions are negatively charged compared to Ga"^^, 
in the very dilute limit they bind these holes forming 
an acceptor level with a binding energy Eb ~ 112meVii& 
As the Mn concentration increases, these acceptor states 
start to overlap and form an impurity band, which for 
even larger Mn concentrations merges with the valence 
band. Though the actual concentration at which the im- 
purity band disappears is not known, according to op- 
tical conductivity measurementSfSiifl and ellipsometry^^ 
this impurity band seems to persist up to nominal Mn 
concentrations as high as Xnom ~ 0.05. Angle resolved 
photoemission (ARPES) data^iSiiiii scanning tunneling 
microscope (STM) resultsi^*i& and the fact that even 
"metallic" samples feature a resistivity upturn at low 
temperatureii suggest that for smaller concentrations 



(and maybe even for relatively large nominal concentra- 
tions) one may be able to describe Gai_a;Mna;As in terms 
of an impurity bandiiSiiSiSSiSi While optical conductiv- 
ity results^- ^" and ellipsometry resultaii are suggestive of 
the presence of an impurity band in moderately doped 
samplesjSiifi an interpretation based on band to band 
transitions is also possible^^^ We remark, however, that 
these materials are extremely dirty^^~ihe mean free path 
is estimated to be of the order of the Fermi wavelength- 
and therefore it is not clear if the latter approach is ap- 
propriate. Also, ARPES data indicate that the chemical 
potential of "insulating" samples lies inside the gapfi^ 
contradicting a band theory-based interpretation of the 
optical conductivity data. 

A detailed understanding of an impurity band model 
begins with the knowledge of a single Mn acceptor 
statei^ The physics of the isolated Mn^^ -I- hole sys- 
tem is well understood:^ In the absence of the Mn^^ core 
spin, the ground state of the bound hole at the acceptor 
level is four-fold degenerate and well described in terms 
of a F = 3/2 state. For most purposes, only the fourfold 
degenerate F = 3/2 acceptor levels need be considered 
in the dilute limit even in the presence of the Mn^"*" core 
spin. As evidenced by infrared spectroscopy)^ the effect 
of the 5 = 5/2 Mn core spin on the holes is well-described 
by a simple exchange Hamiltonian:^ 

Hcxch = GS ■ F , (1) 

with G « 5 mcV.* 

The bound hole (acceptor) states within the F = 3/2 
multiplet are not Hydrodgenic due to a significant d-wave 
component of the bound state wavefunctionii2i2Si2Si2i 
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This d-wave character ultimately comes from the spin- 
orbit coupling in GaAs and has recently been confirmed 
in the beautiful STM experiments of Yakunin et. alm^ 
The anistropy of the orbital structure of the wavcfunc- 
tion leads to dircctionally dependent hopping of holes 
between Mn ions, a splitting of the F = 3/2 level degen- 
eracy, and is expected to strongly influence the magnetic 
and transport properties of dilute GaMnAs.^^ Here we 
study these effects in detail. 

One of the main results of this paper is thus the ef- 
fective Hamiltonian describing strongly interacting holes 
hopping from Mn to Mn. The holes arc coupled to the 
Mn spins via the exchange interaction where 

^cff ^ ^cff ^ ^. (2) 

The first part of this Hamiltonian, Hq^ , describes the 
hopping of the holes from Mn to Mn, and the interactions 
of the Mn acceptor site with the Mn core spin, 

+ G ^ 5, • {cl^F^'' c.) . (3) 

To determine the parameters of Q we shall use the 
spherical approximation^ This approach neglects the 
cubic symmetry of the lattice, but approximates the band 
structure rather well around the top of the valence band 
at the Gamma point which is most relevant at the low 
hole concentrations of interest in the present paper. The 
term Hint accounts for the on-site interactions of holes 
with each other, and in the spherical approximation, 

i i 

The operator cj ^ in the above expressions creates a hole 
at the acceptor level \F = 3/2, F^ — v) at site i, Ni = 
J2:.4,uCi,u, Fi = J2f,,^clf,F'"'ciM, and : ... : denotes 
normal ordering. Here is the fiiy element of the spin 
3/2 matrix. The Hubbard interaction strength Un and 
the Hunds rule coupling Up in Eq. Q can be obtained 
by evaluating exchange integrals and we find Un ~ 2600 
K and Up « -51 K. 

The presence of nearby Mn sites has three important 
effects on the F = 3/2 acceptor state at any particular 
Mn site: (i) The Coulomb potential of the neighboring 
Mn^"*" ions induces a random (from the random relative 
positions of the Mn) shift E of the fourfold degenerate 
states, (ii) Because of the large spin-orbit coupling in 
GaAs, the neighboring atoms also generate an anisotropy 
K and split the fourfold degeneracy of the F — 3/2 state 
into two Kramers degenerate doublets, (iii) Finally, the 
presence of the neighboring ions allow these F — 3/2 
spin objects to hop between the Mn sites. However, this 
hopping t does not conserve the spin F because of the 
spin-orbit coupling. 
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FIG. 1: Spin orientation dependence of the ground state 
energy of one hole on a Mn dimer parallel to the z-axis for 
a spatial separation of = 14 A. The Hamiltonian Q is 
diagonalized exactly with the parameters shown in Fig. 2] As 
the spins are rotated by an angle 9 away from the z-axis the 
energy increases and reaches a maximum at 6 — 90° before 
again decreasing. This indicates the magnetic anisotropy is 
easy axis. 



To determine the parameters of (O , we performed vari- 
ational calculations for a dimer of Mn ions taken to lie 
along the z-axis where Fz is good quantum number. Once 
the parameters of the dimer is in hand and the positions 
of all the Mn are known, the parameters of the Hamilto- 
nian Q are obtained by simple spin-3/2 rotations. 

To illustrate the power of the approach, and to better 
understand the physical results we obtain from it, con- 
sider the simplest case of 2 Mn impurities and 1 hole. 
Diagonalizing Eq. (j^jl for different orientations of 5*1 and 
S2 with the parameters given in Fig. ^ shows that the 
magnetization has an easy axis anisotropy (See Fig. 
This easy axis anisotropy immediately leads to frustra- 
tion among non-coUinear Mn positions. 

We study the Hamiltonian = ffg^ -I- Hi^t in de- 
tail using mean field theory when Un = Up = and 
also with a variational approach when Un,Uf 7^ 0. We 
study the interplay of disorder and dircctionally depen- 
dent hopping parameters induced by spin-orbit coupling. 
We calculate the temperature dependence of the mag- 
netization, magnetic anisotropics, the spin distribution 
functions measuring the degree of non-coUinearity among 
the spins, the (impurity band) density of states, and the 
dependence of the localization transition on the various 
parameters of our model. Our main results are the fol- 
lowing: Qualitatively similar to our earlier results in the 
metallic regime, we find that the interplay of disor- 
der and spin-orbit coupling results in (i) magnetization 
curves that exhibit linear behavior over a significant tem- 
perature range and (ii) a broad spin distribution function, 
implying highly non-coUinear magnetic states that result 
from spin-orbit induced magnetic anisotropics (iii) within 
our mean field and variational calculation we find a 
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well developed impurity band separated from the valence 
band for active Mn concentration up to ^active ^ 0.015 
with a localization transition fairly robust to interactions. 

In this paper all Mn concentrations x are the active 
Mn concentrations ^ i.e., x — ^active where active Mn are 
defined to be those Mn that contribute to the ferromag- 
netism of the material. Interstitial defects with a Mn 
sitting next to a substitutional Mn may result in a lo- 
cal singlet formation)22i thereby rendering the two Mn 
magnetically inactive since they do not contribute to the 
fcrromagnetism of the material. Thus, the active Mn 
concentration is typically less than the nominal Mn con- 
centration. 

The interstitial Mn also compensate holesS*^ reduc- 
ing the number of itinerant holes. In this paper we 
use the hole fraction /, to relate the hole to the Mn 
concentration as Nh — f Nmu where is the num- 
ber of holes and Nmu is the number of active Mn. Al- 
though the precise value of / is not known, typically, 
/ = 0.1 — 0.3. We thus include the effects of various com- 
pensating defectS)^^i^°i'^^i^^i'^^i^^ such as interstitial Mn 
and As antisites indirectly through the parameter /. 

The outline of this paper is the following. In Sec.Hllwe 
describe the variational calculation used to obtain an es- 
timate of the bound state acceptor wavefunction around 
a single Mn ion. In Sec. IIIII we use the variationally ob- 
tained wavcfunctions to derive and compute the effective 
parameters of the Hamiltonian, Eqs.Q and Q), which 
we then study in detail in Sec. II VI using mean field and 
variational approaches. Finally, in Sec. we discuss the 
main conclusions of our work. Technical details of our 
calculations and various lengthy analytical expressions 
are relegated to the appendices. 

II. VARIATIONAL CALCULATION OF THE 
BALDERESCHI-LIPARI WAVEFUNCTIONS 

In order to study GaMnAs in the dilute limit, we pro- 
ceed stepwise by first obtaining bound state (acceptor) 
wavcfunctions in the single substitutional Mn impurity 
limit and then using these wavcfunctions to obtain ef- 
fective parameters of two- ion and iV-ion Hamiltonians, 
details of which are given in Sec. Illll 

We start from the spherical HamiltonianSSiSSiffi 

2m V ^-^ / e r 

a,l3 

where the central cell correctior!i22ii& 

Vcc = -Voe-^^^' , (6) 

is used to reproduce the experimentally obtained bind- 
ing energies, and therefore reasonable acceptor wavcfunc- 
tions. This affects the parameters given in Fig. 0| of 
the effective Hamiltonian Q. Here tq is a short dis- 
tance cutofi^ for the central cell correction and Vq its 



size. The primary role of the central cell correction © 
is to take into account atomic interactions in the close 
vicinity of the Mn ion. In Eq. ij^l 7 ~ 7.65 is a mass 
renormalization parameter, m is the free electron mass, 
fi ~ 0.77 is the strength of the spherical spin-orbit cou- 
pling in the j = 3/2 band of GaAsf2& and e « 10 is 
the dielectric constant of GaAs. The spin-orbit term 
in Eq. ([SJ couples the momentum tensor of the holes 
Pa/3 — PaP/3 — 5af3 P^/3 to their quadrupolar momen- 
tum, Jap = (j )/2 - 5a/3 j{j + l)/3. This ef- 
fective Hamiltonian gives a relatively accurate value of 
the hole energy in the vicinity of the top of the valence 
band, but is not very reliable for holes with higher en- 
ergy, since then other states not included in the deriva- 
tion of (0) will be mixed into the acceptor state wave 
functions. The Hamiltonian ||3Jl also does not distinguish 
between different crystalline directions. We will discuss 
the implications of these features and other shortcomings 
of the spherical approximation in the concluding section, 
Sec.lVl 

To proceed with the calculation, we note that Eq. © 
can be made dimensionless by measuring distance in 
units of the effective Bohr radius, OcS — h'^ej/e^rn = 
40A, and taking the corresponding effective Rydberg, 
^cff = e'*m/2?i^e^7 = 15.7meV, as the energy scale. In 
our calculations we have used ro = 2.8A and Vq — 3.0 
eV. These values are very close to the numbers used for 
the central cell corrections in Ref.ls^ and Ref./40]. 

With the central cell correction, we obtain the cor- 
rect binding energy of 112 meVi^ However, due to the 
central cell correction lO, a^s is no longer a measure of 
the spatial extent of the wavefunction as it would be for 
a purely Coulomb potential. Instead, the characteristic 
length scale is ^ 10 A, as can be seen in Fig. [3 

When ^ 7^ in Eq. ©, the ground state of a hole 
bound to an acceptor is no longer a state of zero orbital 
angular momentum, L = 0, since the "spin-orbit" term 
will mix in a d-wave, L = 2, component. The ground 
state wavefunction is therefore no longer Hydrogenic and 
hence not spherically symmetric li^iSSiS^ This feature will 
lead directly to the appearance of spin-dependent hop- 
ping terms in Eq. |(2Jl. 

Within the spherical approximation, the total angular 
momentum F = L+j is a constant of the motion and for 
/i 7^ and the ground state has F = 3/2. The wavefunc- 
tion for the ground state can then be written as a sum 
of an s-wave component /o and a rf-wave component go 

1>F.(r) = /o(r)|L = 0,j = |i^=|F,) 

+ go{r)\L^2,j^^,F^^,F^). (7) 

By acting with the Hamiltonian, Eq. lO, on Eq. {Tj) one 
obtains the following set of differential equations to be 
solved for fo{r) and goir): 
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(8) 



where Vr, 



2ma„ 



-Vcc- In order to solve Eq. (jSJ we 
follow the variational approach of Ref. by expanding 
/o and go as 



N 



Mr) = ^A,/,(r) 

i=l 
N 

Y.^^9^{r) 



9o[r 



(9) 
(10) 



where the Ai and Bi are variational parameters to be 
determined and the fi{r) and gi{r) are normalized but 
not orthogonal basis functions 



Mr) 



^{r) 




(11) 
(12) 



with ai — a. In our computations we have taken 
AT = 21, a = 1 X 10"^ and aw = 5 x 10^ as in Ref. [2^ . 
and we also verified that refining the basis set resulted in 
no further improvement. 

To obtain the ground state wave function, we minimize 
the expectation value of the Hamiltonian on the left hand 
side of Eq. (jS)). Using the coefhcients Ai and Bi as vari- 
ational parameters this involves the solution of a simple 
2N X 2N eigenvalue problem. One must, however, also 
take into account during this calculation that the states 
fi and gi are not orthogonal. 

The non-orthogonality of the basis set can be taken 
into account through the computation of the over- 
lap matrices, = /o°° r^fi{r) fj{r) and Sf^ = 
dr r'^gi (r) gj (r) , and the transformation of the orig- 
inal problem to a corresponding new orthonormal ba- 
sis following rather standard atomic physics procedures. 
The radial functions /o (r) and g^ (r) obtained in this way 
are shown in Fig. |21 Note that while the s-wave compo- 
nent /o dominates at short distances, the d-wave com- 
ponent go becomes appreciable for r > lOA. This d- 
wave component is ultimately responsible for the strong 
anisotropy of the hopping and effective spin-spin interac- 
tion. 

Using the radial wavefunctions plotted in Fig. |21 one 
can compute the expectation value of the local spin den- 
sity, (i(r)), around a Mn impurity. Replacing the Mn 
spin for a moment with a classical spin pointing down- 
ward along the 2-axis, a bound hole on the acceptor level 
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FIG. 2: (Color online.) Radial wavefunctions obtained from 
a variational calculation for n = 0.767, the relevant value 
for GaMnAs. For r > ISA, the typical Mn-Mn distance at 
X — 0.01, goir) ~ foir)- From Eq. lO this means the d-wave 
component of the wavefunction is important for typical Mn 
concentrations at x = 0.01. It is thus expected that the non- 
Hydrogenic nature of the wavefunction will strongly affect the 
magnetic and transport properties of dilute GaMnAs, a result 
supported by our numerical calculations presented in Sec. llVI 



will occupy the state — 3/2, provided that the cou- 
pling between the Mn spin and the hole is antifcrromag- 
netic. The spin direction (polarization) of this bound 
hole around the impurity is shown in Fig. |2| Note that 
the polarization direction depends on distance and can 
change sign. Note also that in the absence of spin-orbit 
coupling, fj, = 0, the spin polarization of the hole would 
be just pointing along the z direction, and display RKKY 
oscillations at larger distances (not shown in the figure). 
Detailed expressions for the acceptor state spin density 
are given in Appendix fXl 



III. COMPUTING THE TWO-ION AND iV-ION 
HAMILTONIAN 

Using the variational wave function obtained in SeclTTl 
Eq. |(7J), we now compute the effective parameters of the 
two- ion hopping Hamiltonian, Eq. (|15|l . which will in turn 
allow us to find the parameters of the A^-ion Hamiltonian, 
Eq. (j2Jl, by using spin-3/2 rotations. 

We assume that we have two impurities separated by 
a distance R. We take the quantization axis, z, to be 



5 



<J(r, e)> for 11=161 



A 



o 




FIG. 3: Polarization of a bound hole in the state |_F — 
3/2, Fz — 3/2) in Gai_xMnxAs around a Mn ion (dark arrow 
pointing downward represents the Mn S=5/2 spin). Only the 
direction of the polarization is indicated. The magnitude falls 
off on a scale ~ 10 A, as indicated by Fig. |5| 



along the line joining the two impurities (ions). Neglect- 
ing again the effect of the core Mn spin (for the time 
being) , the full Hamiltonian within the spherical approx- 
imation can be written as 



2 -ion 



where 
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{p^-tiJ2Jo.pPo./3)+Viir) + V2ir^ , (13) 



a,/3 



V^{r} 



+ VU\r-n\), 



(14) 



with fi and r2 the locations of the two impurities. 

Having computed the single Mn hole states, we car- 
ried out a variational calculation to construct the molec- 
ular orbitals for a pair of Mn ions in the approximation 
where we considered only linear combinations of the sin- 
gle impurity ground state wave functions >i2iii For a pair 
of Mn spins the full SU(2) symmetry of the single im- 
purity model is broken. However, the Hamiltonian H13|) 
still possesses a cylindrical symmetry, corresponding to 
the conservation of Fz- As a consequence, the various 
Fz subspaces decouple, and our task reduces to the con- 
struction and diagonalization of 2 x 2 matrices. Further- 
more, time reversal symmetry implies that the two states 
with Fz — ±1/2 and the two states with Fz — ±3/2 re- 
main degenerate. As a consequence, we find that the 
two four- fold degenerate F = 3/2 acceptor states of the 
two Mn impurities are split into four Kramers degener- 
ate doublets. (Details of this calculation are given in 
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FIG. 4: (Color online.) Parameters of the two-impurity 
Hamiltonian Eq. 1151 obtained from the variational study of 
two Mn ions. The arrow indicates the typical Mn-Mn dis- 
tance, dtyp, for X = 0.01 Mn concentration. 



Appendix IBl) Since for typical Mn distances these or- 
bitals are well separated from the rest of the spectrum, 
we shall be satisfied by providing a description of only 
these eight lowest lying states of the "molecule". This 
can be achieved by using the following effective Hamilto- 



(15) 



+ EiR) + Eo ct , 



1=1,2 ^ 



where R = \ri — r2|, describes the hopping of Fz — v 
holes, K is the splitting of the F = 3/2 manifold of 
states generated by the presence of the other Mn impu- 
rity, and E denotes the energy shift of the acceptor state 
(at one ion due to the presence of the other ion) with 
respect to the binding energy of an isolated acceptor, 
Eq « 112 meV. By time reversal symmetry, the hopping 
parameters satisfy t3/2 = t-3/2 and = i-i/2- All 
parameters depend only on the distance R between the 
two Mn sites (see Fig. 0J. The most obvious effect of 
the spin-orbit coupling is that the hoppings 1^/2 and ti/2 
substantially differ from each other; holes that have their 
spin aligned with the Mn-Mn bond are more mobile. As 
we mentioned in the introduction, this leads to an easy 
axis magnetic anisotropy in the effective spin-spin inter- 
actions and to non-collinear magnetism. As indicated by 
the arrow in Fig. 0] at the typical Mn-Mn distance for 
X — 0.01, K and ti/2 can be entirely neglected compared 
to E and t^/2- Therefore, in many cases it is enough to 
keep only the latter two terms in the effective Hamilto- 
nian. 

So far, we have neglected the interaction between the 
core Mn spins S and the acceptor state. It is known from 
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experiments,- that the spectrum of an isolated Mn im- 
purity can be very well described by a simple exchange 
Hamiltonian, -ffcxch — G S ■ F. Furthermore, the separa- 
tion ~ 100 meV of the acceptor state from other excited 
states is much smaller than the experimentally found ex- 
change coupling G w 5 meV. We can therefore safely 
treat the exchange field of the Mn spin as a perturba- 
tion. We remark at this point that the Mn ions are, to 
a very good approximation, in a d"^ state, and valence 
fluctuations on the d-levels seem to be rather small, as 
evidenced by an experimentally observed ^-factor close 
to 2 In this spirit, we take into account the effect of Mn 
core spins through the following simple term, 

He.ch = G E E ■ (4^^'^" C.,.) . (16) 
z— 1,2 

Note that in this expression we neglected interactions 
between the core spins and the hole spin on a neigh- 
boring Mn acceptor level. This approximation is cer- 
tainly justified in the extreme dilute limit, and the above 
Hamiltonian does give a reasonable value for the Curie 
temperature at the concentrations we consider. However, 
additional terms may be important for a quantitative de- 
scription of GaMnAs.^'* 

Finally, let us discuss the hole-hole interaction term, 
Eq. I0J). Again, the on site interaction can be greatly 
simplified due to the presence of SU(2) symmetry within 
the spherical approximation. Since holes are fermions, 
two holes can be placed to the four lowest lying acceptor 
states in six different ways. These six states correspond 
to a fivefold degenerate total spin F = Fi + F2 = 2 two- 
hole state and an F = singlet state. The interaction 
term can be thus written as 

Hint ^UdPd + Us Ps , (17) 

where we introduced the four Fermion operators Pd and 
Ps that project to the F = 2 and F = two-hole sub- 
spaces, respectively. With a little algebra we can rewrite 
these expressions in the form Eq. Q), and we can ex- 
press the Hubbard interaction U m and the Hund's rule 
coupling Up in terms of simple Coulomb integrals (see 
Appendix [Cl for details). 

In the more general case, with three or more impurities, 
we need to know how to generalize the Hamiltonian l|15|) 
to the situation where the impurities do not lie along the 
z-axis. We can derive the parameters of Eq. from the 
results of Appendix iBlbv applying appropriate rotations. 

This can be achieved as follows. Assume that we 
have two Mn impurities at positions fi and r2- It is 
trivial to write the hopping part of the Hamiltonian if 
we quantize the spin of the holes along the unit vec- 
tor n = {sm{9) cos(^), sin(0) sin((/)), cos(6')) connecting fi 
and r2- Denoting the eigenvalues oi F ■ n hy i>, we can 
write the hopping part of the Hamiltonian in the simple 
form 

<°n-M„ = E*'>(^)(<^'^2,.+h.C.) , (18) 



where cj - creates a hole at site i with F ■ fi = v, and R 
denotes the separation between the two ions. We need 
to re-express this Hamiltonian in terms of operators that 
create holes with F quantized along the z-axis. This can 
be simply achieved by noticing that these two sets of 
operators are related by a unitary transformation: 

4 = Y.^U{ri%., cl , (19) 

where U is just the usual spin 3/2 rotation matrix: 

U{n)^e">'^'e''^''y . (20) 

Making use of this transformation we can rewrite the 
hopping term in this standard basis as 

^mT-M„ = E {tl2ci.C2y + h.C.) , (21) 

where the hopping matrix is simply given by 

fi^ ^Y.P{ri)]^-t,{R)[U\ri)]^^,. (22) 

It is much simpler to generalize the spin splitting term 
^ K which can trivially be written as 

^^Z1m„ = E K{R)cl^{{ri-F)l,, - h,,.\c,y . 

Finally, the energy shift term is manifestly invariant with 
respect to the spin quantization axis, 

<n-Mn = E{R) E E 4.^^- ■ (23) 
i=l,2 V 

For a finite number of ions the above perturbations add 
up in a tight binding approach, leading to the effective 
Hamiltonian Q with 

i^r = Ew,)[(^.,-i^)'-5],., (24) 

^ = Y\^in,,)\,,U{P^,){V\n,,%,, , (25) 

and 

-So + E^(^y") ■ (26) 

We remark here that for large distances EiK) scales as 
XjR and therefore, strictly speaking, the latter sum is 
not convergent. This unphysical result of our approach, 
which does not take into account screening, can be reme- 
died in our calculation by introducing an exponential cut- 
off of the order of the Fermi wavelength in Eq. H26|l . 

This completes the derivation of the parameters of the 
general Hamiltonian (PI, aside from the on-site interac- 
tion described in Appendix lUl 
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IV. MEAN FIELD AND VARIATIONAL STUDY 
OF THE EFFECTIVE HAMILTONIAN 

In this section we study the effective Hamiltonian Q in 
a mean field theory'*^ and within a variational calculation 
when the interaction Q is also included i^Sii^ Through- 
out this section we shall treat the Mn core spins as clas- 
sical variables. Our main goal is to study the interplay 
of disorder in the Mn positions and spin-orbit coupling 
of the GaAs host on the magnetic properties of dilute 
GaMnAs. Due to spin-orbit effects in the GaAs host, 
the effective Mn spin-spin interactions are expected to 
be anisotropicjSSiii and these anisotropics are expected 
to be greater for smaller concentrations of Mn ions and 
holesii^ 



A. Computational methods 

Most of our calculations have been performed in the 
absence of the interaction term, Hmt, where wc used 
a simple mean field treatment of the spins,2£ In this 
approximation one has to solve a set of equations self- 
consistently. 

The first one of these equations just expresses the fact 
that polarization of the impurity spin Si is generated by 
the effective field G{Fi) generated in turn by the polar- 
ization of the hole spins: 

{S,) = S^{cothia,)~-) , a, = %S{F,). (27) 
a-i cti' 1 

Te second equation gives the effective Hamiltonian of the 
holes that must be used to compute the thcrmodynamical 
average (F^), 

i 

+ GY,{S,) ■ {clj^'' c,,.) . (28) 

Here the last term simply expresses that a non-zero av- 
erage of {Si) acts as a local field on the holes and tries 
to polarize them. Note that the latter Hamiltonian is 
quadratic. Therefore, once it is diagonalized and its 
eigenf unctions are constructed, we can construct the cor- 
responding density matrix and compute the finite tem- 
perature expectation values {Fi) in a relatively straight- 
forward way, and thus solve the above equations itera- 
tively. 

Although the Hubbard couphng U = Un in Eq. Q 
is rather large, at small hole fractions two holes overlap 
with a small probability, and therefore this interaction 
term is not expected to play a crucial role.'^— To verify 
these expectations, we carried out calculations for the 



interacting Hamiltonian with Un 7^ at T = temper- 
ature. The Hund's rule coupling Up being rather small, 
we neglected this interaction term throughout these com- 
putations. 

A full Hartree-Fock treatment oi U = Un is cumber- 
some: it requires the self-consistent determination 18 ef- 
fective fields at each site, and we typically experienced 
serious convergence problems while trying to determine 
these fields. However, the essential effects of the interac- 
tion term Q can be captured by a simpler approach that 
retains the variational character of Hartree-Fock theory. 
In such a variational approach, we replace the interacting 
Hamiltonian H{Si) by a non-interacting Hamiltonian 

H,ar{^l^, K Si) = Hf {{Si}) - ^ ^1,N, +Y,h, ■ F, , 

i i 

(29) 

where the variational parameters hi and /^i are numeri- 
cally determined by minimizing (for fixed {5*^}) the ex- 
pectation value of the full Hamiltonian (0var|^^|0var)t)ar, 
Eqs. Q and Q), using the ground state 0var of H^ar- 

A T = minimization with respect to the spins Si 
leads to the condition that the spins must be aligned 
anti-parallel to the expectation values of the correspond- 
ing Fi in this variational ground state. Therefore, af- 
ter finding the expectation values {Fi) in the variational 
ground state for a given spin configuration {Si}, we gen- 
erate a new spin configuration by aligning all spins anti- 
parallel to the {Fi)^s. This procedure is then iterated 
with the new values of Si to obtain a self-consistent vari- 
ational solution that includes the effect of interactions. 
In practice, even this restricted approximation is very 
time-consuming because the minimization of the varia- 
tional energy at fixed Si is computationally expensive. 
The procedure outlined above could therefore be carried 
out for only very small sample sizes. Below, we therefore 
present results obtained through a restricted variational 
approach that only uses the variational parameters fii at 
each site. For satisfactory convergence of the variational 
energy minimization step, we slowly crank up Un from 
to its final value in steps of AU — 10 K. 

In our calculations we considered samples of fixed size 
L = 10 aiat and L = 13 aiat where oiat is the length of 
the edge of the FCC unit cell. The effective Hamiltonian 
Q and Q is only expected to be valid in the very di- 
lute limit of Gai_a;Mna;As, so we considered only active 
Mn concentrations x = 0.005, 0.01 and 0.015. The valid- 
ity of our approach can be checked post-facto by noting 
that the high energy tail of the impurity band has fairly 
small overlap with the valence band density of states for 
these concentrations, as seen in Figs. I8I10I Compensa- 
tion effects have been taken into account through the 
hole fraction parameter, /. Although this parameter is 
not precisely known for low-concentration samples, we 
used the values / = 0.1 — 0.3, typically assumed in the 
literature. 

In order to control the amount of disorder, we intro- 



8 



duced a screened Coulomb repulsion between the Mn 
ions and let them relax using zero temperature Monte 
Carlo (MC) simulations as described in Ref. [sD]- We 
found that the Mn ions relax to their long time config- 
uration approximately exponentially fast with a charac- 



teristic relaxation time t 



clax 
MC 



2.5, and that for long 



times the Mn ions form a regular BCC lattice with some 
point defects. Such calculations are not meant to model 
real defect correlationsS^iSS in GaMnAs, but rather to 
help understand how the disorder in the material affects 
its physical properties, especially when random ion po- 
sitions are important as they are for small x and small 
carrier concentrations 

Once the Mn positions are fixed in a given instance, 
the mean field equations derived from Q are solved self- 
consistently:"*^ We usually start the iterative procedure 
from a configuration where all Mn spins are aligned in one 
direction. We used periodic boundary conditions and im- 
plemented a short distance cutoff in the hopping param- 
eters of Eq. © which corresponds to about 8 neighbors 
for each Mn. The use of this cutoff is justified by the 
observation that our molecular orbital calculations are 
only appropriate for "nearest neighbor" ion pairs, and in 
reality, holes can not hop directly over the first "shell" of 
ions. 




FIG. 5: (Color online.) Top: The dependence of the spin 
distribution function, P{cos{8)), on the on-site interaction, 
U — Un, and the Monte Carlo time for L = lOaiat, x — 0.01, 
and / = 0.30. We averaged over 10 samples. Even with in- 
teractions and at large Monte Carlo times (small disorder) 
the spin distribution function remains broad. This is consis- 
tent with the strong reduction of the saturation magnetization 
(~ 60%) observed in our calculations, independent of Monte 
Carlo time. Bottom: Dependence of the spin distribution on 
the hole fraction f,U = 0, obtained after averaging over 100 
samples. 
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FIG. 6: (Color online.) Magnetization as a function of tem- 
perature, Mn concentration x, and Monte Carlo time tMC 
for different hole fractions /. Here L — lOaiat, (7 = 0, and 
100 samples are averaged over. Top: Hole fraction / — 0.1. 
Bottom: Hole fraction / = 0.3. In both cases, as the Monte 
Carlo time increases, for fixed x, the saturation magnetization 
at zero temperature increases. For both values of / the curves 
remain linear over a fairly wide temperature range, much the 
same as for experimentally measured curves for unannealed 
GaMnAs. The saturation magnetization never reaches more 
than ~ 60% of the fully saturated value. This is consistent 
with the wide spin distribution shown in Fig. |^ and indicates 
that the ferromagnetism is non-coUinear. 



B. Results 

1. Magnetization 

Similar to the metallic case within the spherical 
approximatioufSS we find a ferromagnetic state with a 
largely reduced magnetization, ~ 0.4 for L — 

lOaiat. (See Figs. El and □) We find that this re- 
duction is largely due to spin-orbit coupling, and that 
cos{9i) — di ■ n, (where n is the direction of the ground 
state magnetization vector) has a broad distribution, 
P(cos(6')), quantitatively similar to earlier results ob- 
tained in the metallic case using the 4-band spherical 
approximation in the completely disordered case.^^ The 
interaction Hamiltonian appears to have a negligible 
effect on the spin distribution. Also, relaxing the Mn im- 
purities to form a regular BCC lattice as described above 
appears to have little impact on the spin distribution. We 
checked that this result is valid at least for / = 0.1 — 0.3. 
This is qualitatively different from the metallic case which 
showed a significant sharpening of the distribution func- 
tion as the Mn positions became more ordered, and a 
corresponding increase of the saturation magnetization 
to an almost fully polarized state»^ 

The magnetization for [/ = is shown in Figs. and 
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FIG. 7: (Color online.) Magnetization as a function of tem- 
perature, hole fraction /, and Monte Carlo time tMC for dif- 
ferent Mn concentrations x. Here L = lOaiat, U = and 100 
samples are averaged over. Compare with Fig. |n] Top; Mn 
concentration x = 0.005. Middle: Mn concentration x = 0.01. 
Bottom: Mn concentration x = 0.015. The general trend is 
the same as in Fig.|S| Longer Monte Carlo times lead to larger 
zero temperature magnetizations. The saturation magnetiza- 
tion is roughly independent of Mn concentration x. 



13 The curves indicate that the system never reaches the 
fully polarized state, even for long Monte Carlo times. 
However, as the disorder is reduced the saturation mag- 
netization increases from « 20 — 25% to 40 — 50%. The 
magnetization curves exhibit linear behavior over a large 
temperature range, qualitatively similar to experiments 
on disordered samples. 

Unfortunately, since the numerical calculations are 
rather demanding, we could not perform a proper finite 
size scaling analysis. Therefore, while our calculations 
suggest that the ground state of our model is ferromag- 
netic, we cannot exclude the possibility of a paramagnetic 
or spin glass state for these small concentrations. 



2. Density of States 

We compute the DOS from the Hamiltonian (O and 
in the interacting case Hyar- The results are shown in 
Figs. ISITUI 

Fig. |S1 shows the dependence of the DOS on doping 
X for fixed MC time and Un = Up = 0. The total 
number of states is proportional to x. The overall shape 
is fairly independent of x, over the range of x consid- 
ered here, which shows a peak near the binding energy, 
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FIG. 8: (Color online.) Top: The dependence of the density 
of states on doping x for L = 10 aiat, / ~ 0.50, and tMC ~ 0. 
Data is the average of 50 sample realizations. The DOS is 
normalized to the volume of a unit cell, so the total number 
of states is proportional to x. The half-width of the impurity 
band ranges from 1000 — 2500 K and is centered around -1100 
K, the binding energy of a hole at an isolated Mn. The shape 
of the density of states changes little with the Mn concen- 
tration, X, over the range of values shown. The value of the 
Fermi energy is « —5000 K. For comparison, the valence band 
density of states is also shown. Bottom: The dependence of 
the participation ratio on doping x. Data is averaged over 
50 samples. Larger samples have larger values of the PR for 
delocalized states, while for localized states the PR is inde- 
pendent of system size for fixed x, tMC = 0. The energy value 
that separates L-dependent PRs from L-independent PRs is 
the mobility edge. This depends on x and is larger for larger 
X. For the disordered samples here, the mobility edge is not 
too sharp and lies in the tail of the density of states. 



El, « —HOOK, of the isolated Mn-|-hole system and a 
half-width of 0.1-0.25 eV. The impurity band slightly 
overlaps the valence band DOS. However, comparison 
with the valence hole density of states suggests that at 
concentrations x ^ 0.015 a well-formed impurity band 
may still be present, and it might persist to higher con- 
centrations. Indeed, this scenario seems to be supported 
by many experimentsi ?'^?i^^'^^i^?'^'*i^^'^? 

Fig. ini shows the dependence of the DOS on MC time 
for fixed x. For tMC = the Mn ions are completely 
random while for Imc = 10 the Mn ions form a nearly 
perfect BCC lattice with a few point defects. The main 
effect of disorder, mostly due to the random Coulomb 
shift of Ei in Eq. H26|l. is to broaden the impurity band 
DOS. In the ordered case, the width of the impurity band 
is determined by the value of the dominant hopping pa- 
rameter, ^t the typical Mn separations. 

Fig. EH shows the effects of the interactions on the 
DOS. Within the variational calculation, the absolute 
scale of the quasi-particle energies is not given. How- 
ever, as shown in Fig. ^1 the overall shape of the single 
particle density of states and the energy-dependent par- 
ticipation ratio are almost identical to what we found in 
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FIG. 9: (Color online.) Top: The dependence of the density 
of states on Monte Carlo time for U = 0, x — 0.01, / = 0.30, 
and L — lOaiat. Data is the average of 50 sample realiza- 
tions. MC time tMC = means that the Mn positions are 
completely random; for Imc = 10 the Mn ions form a nearly 
BCC lattice with a few point defects. The main effect of dis- 
order is thus to broaden the impurity band. The width of the 
impurity band is proportional to the value of the dominant 
hopping parameter, ^3/2, at typical Mn separations as can be 
seen from Fig. 01 Bottom: The dependence of the PR on the 
MC time for x = 0.01 and / — 0.30. Data is averaged over 
50 samples. The mobility edge also moves up to higher en- 
ergy values for more ordered Mn configurations and nearly all 
states become delocalized. 



our calculations performed for the non-interacting model. 

In order to gain information on transport properties 
of the holes, we turn to an analysis of another quantity, 
the participation ratio, from which finite size scaling will 
be able to tell us which states of the impurity band are 
localized and which states are delocalized in the impurity 
band. 



3. Participation Ratios 

The participation ratio, PR = E,:(I]q 
measures the degree to which wavefunctions are local- 
ized. If states are completely delocalized, the single par- 
ticle wavefunction ipia will be spread equally over all sites 
making the PR system-size dependent because the wave- 
function must be normalized to unity. Thus, the PR 
grows with system size for delocalized states while the it 
remains 0(1) in the thermodynamic limit for localized 
states. 

Fig. |S1 shows the dependence of the PR on x and sys- 
tem size L for fixed disorder. Larger samples have larger 
values of the PR for delocalized states while for local- 
ized states the PR is L independent. The energy joining 
the two regimes is the mobility edge. It is impossible to 
determine the precise position of the mobility edge from 
our numerics, but in both cases, the Fermi energy appar- 



FIG. 10: (Color online.) Top: The dependence of the density 
of states on the on-site interactions, U — 2600 K, for x = 0.01, 
/ = 0.30, and tnjc- The effect of interactions on the DOS is 
minimal: The overall shape remains largely unchanged by the 
interactions. Bottom: The dependence of the PR on the on- 
site interactions, U, for x — 0.01, / = 0.30, L — lOaiat, and 
tMC ~ 0. Data is averaged over 10 samples. The behavior of 
the participation ratio follows roughly that of the density of 
states shown in Top; There is little shape change compared 
to the non-interacting case. 

ently lies in the region of delocalized states, indicating a 
localization transition in the impurity band itself. 

Fig. 121 shows the dependence of the PR on disorder. 
For small disorder, nearly all states become delocalized 
and similar to the the disordered case the localization 
transition occurs in the impurity band. 

Fig. 1101 shows the dependence of the PR on the on-site 
interactions in Eq. The behavior of the PR ratio 

roughly follows that of the DOS shown in Fig.^1 There 
is little shape change with the interactions, and the re- 
sult looks very similar to the non-interacting case. Thus, 
the relation between the mobility edge and the Fermi 
energy remains essentially unchanged implying that the 
localization transition is robust to reasonable on-site in- 
teractions. 

To summarize the result of this section, we find that 
the chemical potential Hes deep (~ 0.5 — 0.7 eV) inside 
the gap. From the PR data, it appears that the chemical 
potential is in the vicinity of the mobility edge, a regime 
where our model is probably more reliable. This suggests 
that the localization phase transition in Gai_2:Mnj:As 
could happen inside the impurity band and that the fer- 
romagnetic phase for smaller Mn concentrations is gov- 
erned by locahzed hole statesi2i2ai2iiMSiiS4I 

V. CONCLUSIONS 

Starting with a single Mn acceptor state in GaMnAs, 
we derived an effective Hamiltonian for Gai_a:Mn2:As 
valid in the dilute limit, where Gai-ajMn^jAs can be de- 
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scribed in terms of spin F = 3/2 polarons hopping be- 
tween the Mn sites and coupled to the local Mn spins. 
We estimated the parameters of this model from micro- 
scopic calculations using both a variational approach and 
an exact diagonalization for a pair of Mn ions within the 
spherical approximation. Our approach treats the ex- 
tremely large Coulomb interaction in a non-perturbative 
way, and captures the effects of strong spin-orbit cou- 
pling, and disorder. We find that due to the large spin- 
orbit coupling of GaAs, the hopping matrix elements of 
the polarons depend on their spin direction. 

We studied the above effective Hamiltonianon using 
mean field and variational methods, also including the ef- 
fects of interactions between the holes. We find that the 
spin-dependent hopping generates frustration and is ul- 
timately responsible for the formation of a non-coUinear 
magnetic state for small active Mn concentrations. The 
existence of such non-collinear ground states is indeed 
supported by experiments, where a substantial increase 
in the remanent magnetization is found upon the appli- 
cation of a relatively small magnetic field in some unan- 
nealed samples. 

Our calculations also support the existence of an im- 
purity band, and a metal-insulator phase transition in- 
side this impurity band for these small concentrations 
of active Mn ions, in agreement with angle resolved 
photoemission (ARPES) datajiSiiii^ scanning tunneling 
microscope (STM) resultSfi^ii^ and optical conductivity 
measurements. ^'^^ 

The main advantage of our approach is that it provides 
a clear description of the most important physical ingre- 
dients needed to describe dilute Gai_a;Mna;As, while it 
treats the extremely large Coulomb potential of charged 
substitutional Mn ions non-perturbatively. While the re- 
sulting effective Hamiltonian given by Eqs. ^ and Q is 
relatively simple, it captures many of the physical prop- 
erties of Gai_a;Mna;As, and can serve as a starting point 
for field theoretical computations of other physical quan- 
tities of interest such as optical conductivity, spin wave 
relaxation rate, conductivity or (anomalous) Hall resis- 
tance. 

Though the parameters of our effective Hamiltonian 
have been determined from microscopic model calcula- 
tions, they are only approximate: while the spherical ap- 
proximation used is able to reproduce the spectrum of 
a single acceptor rather well, it certainly overestimates 
the effect of spin-orbit coupling and the width of the im- 
purity band. A direct comparison of the parameters in 
Fig. ^ with those obtained from a more accurate six- 
band model calculation shows some important quanti- 
tative differencesJ^^ This comparison reveals that while 
for r > I3A the effective Hamiltonian © is indeed a 
good approximation in form, the hopping parameters are 
smaller by a factor of two compared to the ones obtained 
from the six band model variational calculation. More- 
over, for r < 13A, the six-band model gives ^3/2 ~ ^i/2j 
suggesting that spin anisotropy is much smaller than that 
obtained from the spherical model. Furthermore, for 
shorter Mn separations the effective model turns out to 



be a rather poor approximation^ 

In summary, based on microscopic calculations, we 
constructed a many-body Hamiltonian that is appropri- 
ate for describing Gai_a;Mna;As in the very dilute limit, 
and estimated its parameters. We find that the hopping 
of the carriers is strongly correlated with their spin. This 
spin-dependent hopping is crucial for capturing spin- 
orbit coupling induced random anisotropy terms, or the 
lifetime of the magnon excitations. Our calculations sup- 
port the presence of an impurity band for ^active ^ 0.015 
active Mn concentration. 
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APPENDIX A: EXPRESSIONS FOR ANGULAR 
DEPENDENCE OF INDUCED HOLE 
POLARIZATION 

With the wavefunctions Q in hand, we can calcu- 
late the average hole spin density around an isolated Mn 
impurity, (j(r, 0, (/>)), which reflects the partial d-wave 
character of the Baldereschi-Lipari wavefunctions. As an 
example, consider {jz{'r,0,4>)) f^=3/2- Using the angu- 
lar momentum addition rules we can express the orbital 
parts of the wave functions in Eq. ITJ as 

\L = Q,,^\^F=\,F.^\)^YS\\), (A1) 
\L = 2,j = %F=l,Fz = l)^ (A2) 




where the Yj"* denote the spherical functions, and the 
second terms stand for the spin part of the j = 3/2 wave 
function. Thus the full wave function reads 




(A3) 

Taking the expectation value of j\\ = jz in this state 
gives, along with analogous calculations in the other 
states and for the perpendicular component of the spin, 
j_L = cos(0) -f sin(0) jy, 
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(j||(r))_F, = ±3/2 
(J-LW>_F, = ±3/2 

(j||(r)>_F,=±i/2 

(j±(f^))_F, = ±l/2 



±^ [/oW + /o(r).go(r)(3cos2(0) - 1) + .go(r)^cos4(0)] , 



_^ 3 
47r 



(Mr) 



50 W 



(3cos^(0) - l))go(r) +goWsin^(f?) 



siii(0)cos(0) , 



±- 



1 

'Stt 



Mrf - /o(r)go(r)(3cos2((?) - 1) + ^^^(5 + 12cos(20) ~ 9cos(40)) 



(/oM-^(3cos2W-l))5oW 



sin(6')cos(6l) . 



(A4) 
(A5) 

(A6) 
(A7) 



APPENDIX B: TWO-ION PROBLEM 



while for the Fz ~ ±1/2 subspace we have 



Here we derive the parameters of the effective Hamil- 
tonian H15() using the molecular orbitals for a pair of Mn 
ions.^^ Since the exchange interaction with the Mn core 
is much less than the binding energy of the holes, and 
the on-site interaction energy, we neglect its effect on the 
parameters of the effective Hamiltonian ()15|l . The local 
field created by the Mn core spin on the holes is later 
treated self-consistently in a mean field and variational 
calculation described in Sec. IIVI 

We solve the problem in the eight-dimensional sub- 
space spanned by the F — 3/2 acceptor states centered 
on each impurity obtained through the variational calcu- 
lations of Sec.|nJ As we discussed in the main text, within 
the spherical approximation used throughout this paper, 
Fz is conserved if the two impurities are aligned along 
the z axis. In this case the sectors of different Fz decou- 
ple. Furthermore, because of time reversal symmetry, the 
overlap matrices S (see Sec. ^ and Hamiltonian matrix 
elements are identical for Fz = ±3/2 and for Fz = ±1/2. 
In the Fz = ±3/2 sector these are given by 



1 a+ 



and 



//-(3/2 



Eq + ei 
63 + a+ Eq 



63 + a- 



V Eq 
ei 



(Bl) 



S'(l/2) = 
ij(l/2) ^ 



1 a_ 
a- 1 



Eq +62 64 + a_ Eq 

64 + a_ Eq Ea + 62 



(B3) 



(B4) 



The two columns of these matrices correspond to the two 
Mn sites, and the constants a±, and ei, . . . 64 denote var- 
ious matrix elements between the wavefunction of a hole 
at site 1 and a hole at site 2. The explicit formulas for 
these quantities are given below. Eq is the ground state 
energy of the single bound hole as determined in Sec. 
Using Eq. Q, expanding the angular parts in spherical 
harmonics and then rewriting the expressions in cylindri- 
cal coordinates, we have r = \J + z'^, with p the radial 
coordinate. To simplify our expressions, we introduce the 
notations /o = fo{r(p,z)) (likewise for go), z = z — zq, 
f = p'^ + {z — zo)^i and /o = /o(^) (and likewise for 
5o = 3o(^))i with zo the distance between the two impu- 
(B2) rities, and express the above matrix elements as 



a± = 



ei,2 = 



63,4 



pdp / dz 



if/o±^(3^-l 



/•oo poo 

/ pdp / dz V2{f) 

Jo J -00 

pdp / dz V2{f) 



^4(3j-l 



3 ~ p I ~ - 2 

7; 90 90 [ZZ+-P 



3 2 P / 2,^2 



2^0^ +IP- 



3 . p' f . 1 2 
2^0 90 ^ i^^+^P 



(B5) 
(B6) 
(B7) 



where V2 = Vi (f) is given by Eq. H14|) . It should be kept in mind that a± , the hole binding energy Eq and the four 
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Ci all depend on the spherical spin-orbit strength /i, and 
must be evaluated numerically. These parameters are 
shown in Figs. and El Having these parameters at 
hand, we can simply determine the effective parameters 
^1/2; ^3/2 1 E^J^d E in Eq. (|15|l by equating the spectrum 
of the two Mn ions with that of the effective Hamiltonian 
Eq. H15|) . In this way we obtain 



13/2 



^3 
tl/2 

K 
E 



a+ei - 63 

a\-l 
a_e2 — £4 
- 1 

1 / a+es — ei a_e4 — 62 
2 



4 



al - 1 

1 / Q+e3 - ei a_e4 - 62 



(B8) 
(B9) 
(BIO) 
(Bll) 



These parameters have been plotted in Fig. ^ 
Ir 
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FIG. 11; (Color online.) Wave function overlaps for the two 
states on site 1 and site 2. The overlaps are computed from 
Eq. (B5). 



APPENDIX C: DERIVATION AND 
EVALUATION OF ON-SITE INTERACTIONS Un 

AND Uf 

In the dilute limit it is important to include the effects 
of interactions between holes. Here we only consider the 
on-site interaction of the holes which dominate all other 
interactions due to the localized nature of the molecular 
orbitals. 

In second quantized form the interaction between two 
holes is 
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FIG. 12: (Color online.) Interaction overlaps for the interac- 
tions defined in the text, Eqs. (B6) and (B7). 



where U f^j^j^j^ denotes the usual Coulomb integral 



C^/l,/2./3./4 = E / ^'^Y '^'^"*}l(^^'/^)'^'/4(^^'/^) 



e r — r 



$}^(r',z.)<i>/3(r',z.), (C2) 



and where we have again restricted ourselves to the same 
F = S/2 subspace, and correspondingly the eigenvalues 
of the z-component of F^, fi, may take on the values 
±3/2 and ±1/2. Here n, v are the eigenvalues of jz- The 
wave functions $ (r , have been determined previously 
with the variational calculation outlined in Sec^ (See 
Eq. and Eq. (|A3|) for an illustration of how the angu- 
lar dependence of (r, [i) is obtained. A simple projec- 
tion of (j = 3/2,^2 = M = 1/2| into Eq. (|X3|) picks out 

$3/2(r, 1/2) = -g^(r)^Y^{d,^), for example.) 

Fortunately, we do not have to compute all these ma- 
trix elements if we rewrite Eq. (|C1(I in terms of two- 
hole scattering processes and exploit rotational symme- 
try. Two holes can only take an F = or an F = 2 
configuration within the ground state multiplet because 
of the Pauli principle. One can verify by direct evalua- 
tion that the F ^ = 2 two-hole state is created by the 
following operator from the vacuum, 



''1/2 ''3/2 ' 



(C3) 



while the F = 2 states of lower Fz can be produced 
by applying the lowering operator. The corresponding 
operators read: 



/2 ''3/2 ' 



/l>/2,/3,/4 



V2' - 



3/2 ''3/2 



'--1/2 ''1/2 



(C4) 
(C5) 
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D-2 



''-3/2 ''1/2 ' 
-3/2 ''-1/2 



Likewise for the sole F = operator we get, 



/7T*>'^-l/2 '-I 



V2 



1/2 



„t t N 

''-3/2 ''3/2'' 



(C6) 
(C7) 

(C8) 



Since these operators transform as F = and F = 2 
tensor operators under SU(2) rotations, the interaction 
Hamiltonian must have the form 



J?on-sito = Ud ^ DjjiDm + UsSqSq 



(C9) 



We can, however, use instead of the decomposition above 
the following two SU(2) invariants too: 



2 2 



(CIO) 



where : . . . : denotes normal ordering and N = ''/''/' 

and F = J2f f '^\^fj"^f' denote the number of holes 
and their total spin operator. It is easy to determine the 
relation of the constants Up and Un to Ud and C/s if one 
rewrites Eq. (jClOp using the identities 



F^ : — F^ 
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-N 



(Cll) 



and compares the action of Eq. (|C9p and Eq. ljC10|) on 
the N = 2 and F = 0, 2 states. This simple algebra gives: 



Un = 
Up = 



- Us 
4 

Ud-Us 



(C12) 
(C13) 



By comparing the matrix elements of Eq. I|C9|I to the 
matrix elements of Eq. IClfl . we can evaluate Ud and 
Us in terms of the Uj.j^-j^j^ = Uf.j^j^j^ wh ich in 
turn allow us to evaluate the Um and Up Eq. (|C10|I . 
Carrying out this calculation, one obtains 



Ud 
Us 



= 2 U^ 



- U: 



1.1 

' 2 ' 2 ' 



ii 1- 1 ^ 

2 ' 2 ' 2 ' 2 



Ui 



1 1 
' 2 ' 2 



Ud 



(C14) 



(C15) 



and therefore, 
Un 



3 I t/s 1.3 1 

\ 2 ^ 2 ' 2 ' 2 



Us 1.1 3 
2 ■ 2 ■ 2 ' 2 



Up 



Ui 



-( 



+ Ui 1 



(C16) 



U. 



3 1.3 1 

2 ' 2 ' 2 ' 2 



U3 



- Ui _i + Ui 1 . (C17) 

2'2'2-2 2'2'2'2/ ^ ^ 

To obtain a numerical value of Upf and U p we 
must determine the matrix elements [/3/2,i/2;3/2.i/2i 

^3/2,1/2:1/2,3/2; t^l/2,- l/2;l/2,-l/2 and ?7i/2,-l/2;-l/2,l/2 

by evaluating the integrals in Eq. ljC2|l . These integrals 
depend on the radial wave functions that we evaluated 
variationally in Sec. Inland are material (parameter) spe- 
cific. In order to evaluate the integrals in Eq. (|C2p the 
^fj..fi{r) must be decomposed into spherical harmonics. 
Various products of spherical harmonics appear in the 
integrand. The integrals can be evaluated by making use 
of the important formula 



47r 



EE 

1=0 m=-l 



(-1) 



m 



r> ) 21 + 1 



(C18) 

where ft (ft') is the angle of f (r^). Here r> (r<) is the 
greater (lesser) of r and r'. With this formula, most of 
the integrals vanish and the few remaining integrals yield 



Us J_.3 1 

2 ' 2 ' 2 ' 2 



Ua 1.1 a 

2 ' 2 ' 2 ' 2 



Ui _i 1 _i 

2 ' 2 ' 2 ' 2 

Ul 1 
2 ' 2 ' 2 ' 2 







ih ~ h) , 



2/2 



(/1+/2) , 



(C19) 

(C20) 

(C21) 
(C22) 



where the prefactor gives the energy scale of the interac- 
tion, 



eaoff 



31.6 meV , 



(C23) 



and Ii and I2 denote the following integrals: 

/•OO /"OO -1 

h = / r^dr / r'Hr' — {fl{r)+gl{r)) 
Jo Jo ^> 

xifoir')+9lir')), 



(C24) 



/2 = / r'dr / r'2dr'-^/o(r)go(r)/o(r')ffo(r') . 
/o Jo 

(C25) 



Evaluating these integrals one obtains Un = 2570 K and 
Up = -51 K. 
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